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^ ' Abstract 

I , We update the Standard Model predictions of the anomalous magnetic moment of the 

Q \ muon, a^ = {g — 2)/2, and the value of the QED coupling at the Z-boson mass, incorpo- 

1^ ' rating the new e'^e~ — > vrvr data obtained by CMD-2 and KLOE, as well as the corrected 

SND data, and other improvements. The prediction for = 11659180.4(5.1) x 10^'^'' 



^ ■ is about 3 X 10 ^'^ lower than before, and has a smaller uncertainty, which corresponds 

to a 3.4 cr deviation from th 
a(M|)-i = 128.937 ± 0.030. 



^ ' to a 3.4 cr deviation from the measured value. The prediction for the QED coupling is 



It is important to predict the anomalous magnetic moment of the muon, = {g^ — 2)/2, 
and the value of the QED coupling on the Z pole as precisely as possible, in order to test 
the Standard Model and to probe New Physics. For the first quantity, we note that the 
Brookhaven experiment gives the average of the measurements of the yU"^ and /i~ anomalous 
magnetic moments to be [T] 

a;^^P = 11659208.0(6.3) x 10"^°. (1) 

If a statistically significant deviation, no matter how tiny, can be definitively established be- 
tween the measured value a^^P and the Standard Model prediction, then it will herald the 
existence of new physics beyond the Standard Model. In particular the comparison offers valu- 
able constraints on possible contributions from SUSY particles, see, for example, the reviews 
in [2]. The second quantity, the QED coupling at the Z boson mass, M^, is equally impor- 
tant. The uncertainty in its value is one of the major limiting factors for precision electroweak 



physics. It hmits, for example, the accuracy of the indirect estimate of the Higgs mass in the 
Standard Model. 

The Standard Model (SM) prediction of the muon anomalous magnetic moment, a^, may 
be written as the sum of three terms. 



«r=«r"+ar+«r. (2) 



The QED contribution, which includes all the photonic and leptonic loops, and the EW contri- 
bution, which includes the loops involving the W, Z or Higgs bosons, are known very accurately: 
qQed ^ (116584718.09±0.16) x lO^^^ [3] and af^ = (154±2) x 10"" [4j. The main uncertainty 
lies in the final term, which involves the hadronic loop contributions. This term may, itself, be 
sub-divided into three parts 

had _ had,LO , ^had,NLO , ^had.l-b-1 /q\ 

At present the hadronic vacuum polarisation contributions, ahad,LO,NLO^ cannot be calculated 
sufficiently accurately from first principles [5j, but instead are evaluated using dispersion inte- 
grals over the measured cross sections for e^e~ — > 7* — ^ hadrons. For our 2003 predictions [6] 
we found 

^had,LO ^ (692.4 ± 5.9exp ± 2.4,,d) X 10-10 , (4) 
^had.NLO _ (-9.79 ± 0.09e.p ± 0.03rad) X 10-^°, (5) 
where the last error corresponds to the uncertainty associated with the radiative corrections to 
the cross section data. The final term in ([3]), a^^'^'^~^^, is the hadronic light-by-light contribution. 
In our previous analysis, we took, in units of 10-1°, eitheiEl 

^had,i-b-i _ 8.0 ±4.0 or 13.6 ±2.5. (6) 

The major uncertainty in the Standard Model prediction of the anomalous moment, a^, 
comes from the contribution of the e~^e^ ti^ti^ channel. Indeed, the e~^e^ ti^ti^ cross 
section data, available in 2003, give the dominant contribution to (jlj) of about (500 ±5) x lO"!*^. 
Since then the situation has improved considerably. In particular, new precise measurements 
have recently become available from the CMD-2 detector throughout the centre-of-mass energy 
range 0.37 < ^/s < 1.38 GeV [101 [IB [12]. Moreover, the KLOE collaboration have made 
measurements of the cross section by the radiative return method over the range 0.6 < ^/s < 1 
GeV [13]. Finally measurements at the SND detector over the interval 0.4 < ^/s < 1 GeV have 
become available [H], and subsequently have been corrected [15]. Clearly all these dataE will 



^ The first is a representative value of several earlier determinations (see, for example, the review in |7]), 
whereas the second value (which was used in the note added in proof in |6i ) was obtained in [8] . In this paper 
we take the second value; it is consistent with the upper bound found in 9J, but see also 'Note added in proof. 

^ There is the possibility of obtaining indirect information on e^e^ — > hadrons in the energy range ^/s < rrir, 
via the conserved vector current (CVC) hypothesis, using the precision data for the hadronic decays of r leptons. 
However there is a sizeable discrepancy between the data from the e~^e~ — s- tt+tt^ experiments and that extracted 
from the n^n^i/ decay mode of the r. This suggests that the understanding of the CVC hypothesis may be 
inadequate at the desired level of precision, see e.g. the recent discussion in [16]. Also the fair agreement between 
the KLOE data and the CMD-2 and SND data strongly supports the e^e~ data against the t data [IT. In 
addition, r spectral function data from Belle show a significant discrepancy compared to ALEPH and, to a 
lesser extent, CLEO, see [18]. We therefore do not include the r data in the present (or in our previous [6]) 
analysis. 
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Figure 1: The lower continuous (red) line shows the contribution of the e+e" 
in the energy region 0.32 < ^/s < 1.43 GeV to a^^^'^'^ 
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size parameter 6. The (orange) band corresponds to the uncertainty in the contribution. All 
the new tt+tt" data are included except those from KLOE, see text. The upper dashed (blue) 
curve, and the error corridor given by the dotted lines, are the corresponding results obtained 
with the data as they were available for our previous 2003 analysis |i6j. 



have an impact on the determination of (7 — 2 of the muon. 

Here we repeat the analysis described in detail in Ref. [6]. That is, we evaluate the dispersion 
relation 

,^^.Lo^(^\' r^^mm^ (7) 



where the kernel K[s) is a known function (see, for example, eq. (45) of |6]), and 



Ris) = ^ . (8) 



The subscript on cTjJg^j = a^^^^{e~^e~ — 7* — > hadrons) is to indicate that we take the bare 
cross section with no initial state radiative or vacuum polarization corrections, but with final 
state radiative corrections; and Upt = ina'^/Ss with a = a(0). 

To begin, we repeat the analysis of [B] including all the new data, except those from KLOE. 
The latter data have a significantly different energy dependence, and we shall discuss their effect 
on the analysis later. We combine the 'bare' cross section data from the various experiments 
for a given channel in clusters of a given size 6, as described in [6]. For the crucial vrvr channel 
we show, in Fig. [H how the contribution a^'^'^° varies as a function of the cluster size 6. (For 
this channel, we have fixed the cluster size in the small energy interval from 0.778 to 0.787 
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GeV containing the p — u interference effects, to be 5 = 1 MeV, thus improving the quahty 
and stabihty of the fit considerably.) The improvement due to the inclusion of the new tttt 
data is immediately clear from Fig. [H Choosing the cluster size 6 to be 3.5 MeV, as in our 
previous analysis [6], we see that the contribution a^^^'^^ is reduced by 4.6 x 10~^°, and that the 
uncertainty in the value decreases from 5.0 x 10~^° to 3.2 x 10~^°. Moreover the result is stable 
to the variation of the choice of the cluster size in the range from 3 to 10 MeV. Furthermore 
the minimum of the clustering fit to the data, divided by the number of degrees of freedom, 
Xmin/d-o.f., is, to a good approximation, equal to 1.1 for all choices of 6 in this interval. It 
is informative to trace the origin of the reduction of 4.6 x 10~^° in due to the addition 
of the new vrvr data. It comes about equally from the intervals 0.32 < ^/s < 0.6 GeV and 
1 < < 1.43 GeV, with only a small reduction in coming from the intervening energy 
range. This becomes clear from Fig. [21 where the data are displayed together with two fits for a 
cluster size 6 = 4.0 MeV (recall from Fig.[T]the insensitivity of the integrated result to variation 
of S in the range 3 — 10 MeV), which correspond to our previous and the present analyses. A 
comparison of the two bands shows that the introduction of the CMD-2 (and SND) data gives, 
on average, a smaller cross section in the energy regions above 410 and above 1030 MeV in the 
upper and lower plots, respectively, of Fig. [21 

We now discuss the inclusion of the KLOE vrvr data in the analysis. Note that KLOE 
measures the hadronic cross section via the radiative return method in e^e^ ^ ^ vrvr7 at 
the $ factory DA<I>NE in Frascati. Here the observation of initial state photon radiation at 
various energies allows for a determination of the invariant mass spectrum of the vrvr system. 
This analysis is completely independent of the 'direct scan' measurements of CMD-2 and SND 
at VEPP-2M in Novosibirsk which use a tunable e~^e~ beam energy. Unfortunately, the KLOE 
data [13] have a different energy dependence to the other vrvr data sets, especially when compared 
with the recent CMD-2 and SND analyses, see e.g. the discussions in [151 [I9]. Our clustering 
prescription allows overlapping data sets to adjust by an overall constant within the systematic 
error of each set, but does not allow for an energy dependent renormalizatioiif|. If we would 
include the KLOE data in the fit, then they would be normalized upwards by nearly two 
percent]^ This in turn would lead to a sizeable effect in the fit, see Fig. [3l the (blue) line is the 
mean value of the fit including the KLOE data, whereas the (red) band is the error band for 
the fit excluding KLOE. In the dispersion integral this would lead to an artificial overall shift 
of a^'^ upwards by about 5 x 10~^°. As the Xmin/d-o.f., for the 'clustering' fit to the vrvr data, 
would increase significantly from the value 1.1 found above, the improvement of the error on 
aji^ after infiation by ^Xmm/^-O-i- would be very small. However the KLOE data [T3| more 
than adequately populate the energy interval of their measurement, 0.60 < y/s < 0.97 GeV, 
so that we can compute reliably their integrated contribution to the dispersion relation of ([7]). 
We find 

a;"'^^(KLOE only, 0.6 -0.97 GeV) = (385.7 ± 4.9) x 10~^°. (9) 

We have also studied the effect of an energy dependent renormahzation, using either a hnear or a constant 
plus logarithmic form; neither of these forms were able to improve the fit with the KLOE data significantly. 

^This is because there are fewer data points in the low energy region lying high (compared to CMD-2) than 
lying low in the central peak and high energy tail region, and because KLOE's quoted systematic error is about 
twice as large as CMD-2's. 
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Figure 2: nir data in the low and high energy p-resonance tail regions (upper and lower panels), 
compared to our fits using a clustering size of 5 = 4.0 MeV. The light (orange) band indicates 
the error band of the fit without the recent CMD-2 (2006) and SND data which became available 
after our previous analysis [6] , whereas the (blue) line and darker (red) band correspond to the 
mean value and error band of the new fit including all data as indicated in the two panels. 
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Figure 3: Mean value of the fit of all vrvr data in the p central region (blue line) compared to 
the fit without the KLOE data (red band), together with data from CMD-2, SND and KLOE 
(as indicated). 



This should be compared to that obtained in the complete analysis (in which the KLOE data 
were omitted) 

a;:"'^°(fit w/o KLOE, 0.6 - 0.97 GeV) = (384.3 ± 2.5) x 10"^°, (10) 

where, as is to be expected, the error is less. Moreover, despite their different energy depen- 
dence, the agreement of the KLOE integrated contribution with that of the fit of all other data 
is excellent. Therefore, for our SM prediction of a^, we combine both contributions for this 
energy interval. Doing so we arrive at 

^jn,LO ^ (^3g4 g ^ 2.2) X 10"^°, (11) 

for the interval 0.60 < y/s < 0.97 GeV. The effect of the KLOE data is to slightly increase the 
magnitude and reduce the error. As a result, the vrvr contribution in the interval 0.32 < ^/s < 
1.43 GeV, becomes (498.5 ± 2.9) x 10'^°. 
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In addition to the inclusion of the new vrvr data, we have made other improvements in the 
analysis. One small improvement that we have made is to evaluate the vacuum polarisation 
corrections to the data using our recent determination [2U] of the running of the effective 
couphng a(g^) in the time-like region = s, which was obtained using our clustered data set. 
This is more consistent than using the evaluation by Jegerlehner [21j that we did previously. 
This is a correction to a correction, and so hardly changes the result. 

A more important improvement is that for the subleading exclusive channels we have in- 
cluded new data from CMD-2 [22j (27r+27r-, 7r°7, 777), BaBar [23] (7r+7r-7r°, 27r+27r-, 37r+37r~, 
7i^TT~ , pp); and from BES [21] for the inclusive hadronic cross section. The effects of the 
new data are summarized in Table [H It lists the contributions to a^^'^'^'^ which have changed 
since our previous prediction. As expected, the main change arises from the inclusion of the 
new 7CTT data: —4.32 in units of 10"^''. This is partially compensated by changes in other contri- 
butions, such as -1-0.84 (due mainly to the new BES data) and -1-0.60 (due to the new leptonic 
widths of J/tp and ip' [25]). Including all the hadronic contributions we obtain 



,had,LO 



(689.4 ±4.2exp±1.8rad) x 10 



-10 



and, similarly for NLO, 



a 



had.NLO 



(-9.79 ± 0.08cxp ± 0.03rad) x 10 
which are to be compared with (jl]) and ([5]) respectively. 



-10 



(12) 



(13) 



channel 


range (GeV) 


X 10^° 


HMNT03 [6] 




0.32 - 1.43 


498.46 ±2.87 


502.78 ±5.02 




0.66-1.43 


46.18 ±0.94 


46.43 ± 0.90 




0.6125-1.43 


6.01 ±0.19 


6.16 ±0.32 




0.60-1.03 


4.54 ±0.12 


4.50 ±0.15 


777 


0.69-1.35 


0.72 ±0.03 


0.73 ±0.03 


inclusive 


1.43 - 11.09 


74.80 ± 2.67 


73.96 ±2.68 


J/tp + ^' 




7.90 ±0.16 


7.30 ± 0.43 



Table 1: Contribution of different channels to compared to the numbers as given in [6]. 



Finally, adding all the terms of ([2]) and ([3]), we obtain the updated Standard Model predic- 
tion for the anomalous magnetic moment of the muon: 

af^ = (11659180.4 ± 5.1) x 10"^°. (14) 

That is, the difference 6af,{= a;fP - af^) is = (27.6 ± 8.1) x 10"^°. Eq. dHj) should be 
compared to our previous prediction (HMNT(03b)), given in the note added in proof in [6], of 

af^ = (11659183.5 ± 6.7) x 10~^°. (15) 

The new data, particularly for e"'"e^ tttt in the centre-of-mass energy range 0.37 < y/s < 1.38 
GeV, have decreased the value of by 3.1 x lO^^*' and reduced the error from 6.7 x 10^^° 
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Figure 4: In the lower half of the plot we compare the experimental value of the anomalous 
magnetic moment of the muon, [g — 2)/2 = a^, with the SM prediction calculated in the text. 
We see that there is now a 3.4 a discrepancy between experiment and the SM prediction. In 
the upper part of the plot we show four earlier determinations [6l [26] of that were made 
before the CMD-2 [TOl [IH [12] , KLOE and SND [ll[T5] e+e" nn data became available. 
The light-by-light contribution used in these determinations varies from (8.6 ± 3.5) x 10~^° in 
DEHZ(03) to the more recent estimate [S] of (13.6 ± 2.5) x 10"^° used in HMNT(03b,06), see 



to 5.1 X 10~^°. Both of these effects increase the disagreement with the measured value. The 
situation is shown pictorially in Fig. [H We now have a discrepancy of 3.4 o", which is larger 
than befor^f|§|. Moreover, note that we use the recent, larger, value [8J of the light- by-light 

^ If we do not use the KLOE data at all, then we have 

flhad.LO ^ (ggg 2 ± 4 3^^p ± l.Sexp) X IQ-^", (16) 

in place of (HH). (The change in the hadronic NLO contribution is invisible to this accuracy.) The corresponding 
value for the total SM prediction is 

af^ = (11659180.2 ± 5.3) x 10^^°. (17) 

The deviation Ja^ from the experimental value a°^^ is 

Sa,, = al^P - a^^ = (27.8 ± 8.2) x IQ-'^" , (18) 

which again corresponds to a 3.4 cr discrepancy. 

6 A preliminary DEHZ analysis [37], also including the new data, finds af-^ = (11659180.5 ± 5.6) x 10"^°, 
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contribution, near the upper limit estimated in [9], which reduces the discrepancy. 

The larger discrepancy, a^^P — a^, which arises from the inclusion of the new tttt data, 
is becoming a more significant indication of New Physics beyond the Standard Model. The 
effect of supersymmetry on can be seen, for example, from the reviews in [2]. Finally, it 
is worth noting that the theoretical error is now below the experimental error on a^. With 
further measurements of the low energy hadronic cross sections underway, which will improve 
the accuracy of a^, the case for improving the measurement of is strong. 

tanj3=10, ii>0, Ao=-300 GeV, m^l71.4 GeV 
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Figure 5: The contours in the minimal SUGRA mi/2, rno plane allowed by the difference 
^exp _ ^SM ^ ^27.6 ± 8.1) X 10"^° at 2-a and 3-a. The other parameters are tan/? = 10, 
Aq = —300 GeV, fi > 0, and the SPS la' reference point is shown by the big dot. Also shown 
are the constraints from the Higgs boson and the chargino masses from the direct searches and 
the charged dark matter (f LSP). 

As an example of the impact of the difference a;fP - af^, we show in Fig. O the allowed 
region of the minimal SUGRA parameters mi/2 and rriQ. We choose the other parameters to 
be tan/3 = 10, Aq = —300 GeV, and /i > 0, so that the commonly studied reference point 
SPS la' [28] appears as the big dot on the plane. We also plotted the constraints arising from 
the Higgs boson mass and the chargino mass from the negative results of direct searches at 
LEP [25], and that from the lightest charged SUSY particle (the region labeled by "f LSP"). 
(To make this figure, we used the program "SuSpect" [29] to calculate the mass spectrum of the 
SUSY particles except for the lightest Higgs boson mass, for which we used "FeynHiggs" [30].) 

The new data also improve the predicted value of the QED coupling at the Z boson mass. 
The value of a(M|) is obtained from [31] 

a-^ = a{0)-^ = 137.035999710(96) (19) 
for which they quote a 3.3cr discrepancy. 
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using the relation 

a{MlY^ = (l - A«iep(M|) - A«had(M|)) (20) 
where Aaiep (M|) = 0.03149769, and Aahad is evaluated from the dispersion relation 

Aa,.d(M|) = -i^P r d. . (21) 

The hare cross section ^"^^^^(s) is defined below (jS]). In comparison with ([7j), dispersion relation 
(I2TI) is less sensitive to the values of Ch^dl-^) the low energy region. We therefore do not 
expect the error on a(M|) to reduce quite as much as that on a|^^*^'^°, when we use the 
improved e+e" — > vrvr data. 

It is conventional to separate out the top-quark contribution and to write 

Aah,d = A45d + Aa*°P, (22) 

where, using rut = 171.4±2.1 GeV, perturbative QCD determines Aa*°P (M|) = -0.000073(02). 
We proceed as in Ref. [6]. Using the new clustered data to evaluate the dispersion relation (pTj) . 
we find 

Aajf2i(^l) = 0.02768 ± 0.00017exp ± 0.00013,ad (23) 
= 0.02768 ± 0.00022. (24) 

This, in turn, gives 

a(M|)~^ = 128.937 ± 0.030, (25) 

which should be compared to our previous estimate of a(M|)~^ = 128.954 ± 0.031. The accu- 
racy is now 23 x 10~^, which is still the least accurately determined of the three fundamental 
parameters of the electroweak theory; AG^/G^ = 0.9 x 10"^ and AMz/Mz = 2.3 x 10~^, where 
Gn is the Fermi constant. 

Note added in proof: After submission of this letter a new assessment of the status of the 
light-by-light contributions has appeared [52]. If we would use their estimate of a^'^'^'^'^'^ = (11± 
4) X 10~^°, the difference a^^^ — (^^ would slightly widen to (30.2 ±8.7) x 10"^°, corresponding 
to a 3.5(7 discrepancy. 
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